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Objective 


Evaluate two Reynolds-stress turbulence models (RSM's) available in the 
FUN3D unstructured CFD code: the SSG/LRR RSM and the Wilcox RSM. 
This work supports NASA’s Revolutionary Computational Aerosciences 
(RCA) Technical Challenge: 


Identify and down-select critical turbulence, transition, and numerical 
method technologies for 40% reduction in predictive error against 
standard test cases for turbulent separated flows, evolution of free shear 
flows and shock-boundary layer interactions on state-of-the-art high 
performance computing hardware. 


Overview 


The FUN3D code 


The turbulence models 


Test cases — simple yet contain relevant flow physics 
— Transonic diffuser 
— Supersonic axisymmetric compression corner 
— Compressible planar shear layer 


— Subsonic axisymmetric jet 


Summary and conclusions 


The FUN3D Code 


e General purpose flow solver and design tool 


¢ Developed by NASA Langley 

¢ Wide variety of numerical schemes, gas models, turbulence models and 
boundary conditions 

e Unstructured grids 

¢ 2°¢_order finite volume, node-centered 


¢ Roe scheme (default) 
— Other methods available 


¢ SA, SST-V, SSG/LRR RSM and Wilcox RSM used 
¢ fun3d.larc.nasa.gov 


me Wind-US Code 
General purpose flow solver 

¢ Developed and supported by NASA Glenn, the Arnold Engineering Development Center 
(AEDC), The Boeing Co. 

e Structured and unstructured grids 

¢ 2°¢_order accurate finite volume, node-centered, Roe (structured) and 
HLLE(unstructured) — default 

e SA, SST-V, EASM models used 


° www.egrc.nasa.gov/winddocs 


Turbulence Models 


turbmodels.larc.nasa.gov 


Spalart-Allmaras (SA) one-equation model 
¢ Standard incompressible version 
¢ No trip term 
¢ p/v=5 freestream boundary condition 


Menter’s shear-stress transport (SST-V) two-equation model 
¢ Vorticity-based production term 


Two-equation explicit algebraic Reynolds stress model (EASM) (shear 
layer case) 
¢ Derived from reduced form of Reynolds stress transport equations 
¢ Similar to the Boussinesg approximation but includes terms that are 
nonlinear in the strain and rotation rate tensors 


Seven-Equation Omega-Based Full Reynolds Stress Turbulence Models 
¢ Wilcox Stress-Omega Full Reynolds Stress Model (Wilcox RSM) 
¢ SSG/LRR-Omega Full Reynolds Stress Model (SSG/LRR RSM) 


Turbulence Models, cont’d 
Seven-equation omega-based full Reynolds Stress models 


SSG/LRR-Omega Full Reynolds Stress Model: 


def VI al 
Tij = —U; Uj 


6 Reynold’s Stress Equations and 1 Length Scale Equation: 


O(pT;; ) ) Faas Up) _ 
Ot Ox;, 


0(pw) d(pwii,) W PPrr _ — Ow p dk dw 
o,. — = a ~ maxiti— — ,0 
Ot ™ Ox, “ok 2 Bu Po v7 ae OX, + oo Ox, OX, 


—pP — ply + pe; — pDy — pMi 


- Blended Speziale-Sarkar-Gatski/Launder-Reece-Rodi pressure-strain model 


1 1 a, 
i =- (G2 + 5 CP ex ) dy a3 Cre (a Aix Ay — 3 Au Ay 5, | + (C3 — C3f Gig Guy JES; 
+ Cyk (tue Sp + Gj Six — 3 Ha Sx 5, + Csk(Gin Wi, + Gy, Wi ) 


Wilcox Stress Omega Full Reynolds Stress Model: 


- Uses a Launder-Rodi-Reece pressure-strain model 


Overview 


The FUN3D and Wind-US codes 


The SSG/LRR and Wilcox Full Reynolds stress models 


Test cases — simple yet contain relevant flow physics 
— Transonic diffuser 
— Supersonic axisymmetric compression corner 
— Compressible planar shear layer 


— Subsonic axisymmetric jet 


Summary and conclusions 


Test Cases 
Transonic Diffuser — Strong Shock Case 
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Test Cases 
Sajben Diffuser — Strong Shock Case 


FUN3D RSM Results 
Wall Pressure 
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Test Cases 
Sajben Diffuser — Strong Shock Case 


FUN3D RSM Results 
Velocity 


y/H, H = local duct height 


x/H = 2.9 x/H = 4.6 
1.2 1.2 
[oe] greets : 
20 
08 2 
oO Experiment o 
064+  —aasase SA 3 
— — SST-V re 
sia — Wilcox RSM = 
— - -SSG/LRR RSM = 
= 
~~ 
T 1 
100 50 Pe 0 50 100 150 200 250 300 1 300 
: u, m/s 
x/H=7.5 
12 > 
- - 
=) =) 
au ou 
= = 
az} as} 
s s 
cs cs 
ws we 
8 8 
T Wl 
= = 
= = 
3 ; 3S T T 1 
-5 300 -50 0 50 100 150 200 250 300 
-0.2 


10 


local duct height 


y/H,H 


local duct height 


y/H,H 


Test Cases 
Sajben Diffuser — Strong Shock Case 


FUN3D RSM Results 
Axial Turbulence Intensity 
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Test Cases 
Sajben Diffuser — Strong Shock Case 


Summary of Results 


¢ The two stress-omega models give very similar results. 


¢ Axial turbulence intensity profiles show better agreement with 
experiment than the SA and SST models. 


¢ The velocity profiles show that the SA model does the best job of 
predicting the separation, however the stress-omega models are better 
at predicting the velocity profiles in the downstream portion of the duct. 


Test Cases 
30° Axisymmetric Compression Corner 
Experiment 


e J. Brown et al, NASA Ames 
¢ Mach 2.85, Re=16x 10°/m 
* Data 

— LDV 


¢ Mean velocities 


|}+——__——-1m 


TEST SECTION 1 


¢ Reynolds stresses WINDOWS 


CUSPED NOSE 


— Surface static pressures OBJECT BEAM 


APERTURE 


Interferometry 
— Schlieren 
Oil flow 


- Dunagan, S.E., Brown, J.L. and Miles, J.B. ,” Interferometric Data for a Shock/Wave Boundary-Layer 
Interaction,” NASA TM 88227, Sept. 1986. 


- Brown, J.D., Brown, J.L. and Kussoy, M.I., “A Documentations of Two- and Three-Dimensional 
Shock-Separated Turbulent Boundary Layers,” NASA TM 101008, July, 1988. 


- *Settles, G.S., and Dodson, L.J., “Hypersonic Shock/Boundary-Layer Interaction Database NASA CR 
177577, April 1991 


- Wideman, J., Brown, J., Miles, J., and Ozcan, O., “Surface Documentation of a 3-D Supersonic 
Shock-Wave/Boundary-Layer Interaction,” NASA TM 108824, 1994 
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Test Cases 
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Mach Contours 
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Xx (cm) 


Test Cases 
30° Axisymmetric Compression Corner 
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Test Cases 
30° Axisymmetric Compression Corner 
Velocity Profiles — Upstream of Flare 
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Test Cases 
30° Axisymmetric Compression Corner 
Velocity Profiles - Downstream of Flare Corner 
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Test Cases 
30° Axisymmetric Compression Corner 
Turbulent Shear Stress — Upstream of Flare 
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Test Cases 
30° Axisymmetric Compression Corner 
Turbulent Shear Stress — Downstream of Flare Corner 


x = 0.433 cm x = 1.732 cm 
2.5 . 
6 Experiment 
SA 
SST-V 
2 SSG/LRR 
Wilcox RSM 


Por 5 -0.01 -0.005 0 0.005 


; 8015 -0.01 -0.005 0 0.005 
<Uuv>/U) 


<uv>/u> 


x = 3.464 cm x = 5.896 cm 


9006 -0.004 -0.002 , 0.002 0.004 19 
<UV>/U 


o 
$006 -0.004 -0.002 5 0 0.002 0.004 
<uv>/U,, 


ny, 


8 


Test Cases 
30° Axisymmetric Compression Corner 


Summary of Results 


The Wilcox and SSG-LRR RSMs behaved quite differently. 
The Wilcox RSM and the SST-V model have similar behavior 


The Wilcox RSM predicted the correct pressure rise on the compression 
surface, whereas the SSG-LRR RSM significantly under-predicted the 
pressure rise. 


The SA model did the best job of predicting the separation location and the 
pressure rise. It also did the best job at predicting the velocity profiles. 


The Wilcox RSM may have an advantage at predicting the shear stress 
profiles. 


Conclusion 
While the Wilcox RSM may offer some slight benefits in predicting the shear stress 
profiles for this case. The SA model gave the best results overall. The SSG-LRR 
RSM performed poorly. 
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Test Cases 
Compressible Mixing Layer 


Experiment 


¢* Goebel, Dutton, & Gruber- Univ.of Illinois (1991) 
* Test Case 2, Convective Mach No., M, = 0.46, Re = 12x10°/m 
¢ Data available: 

¢ LDV Mean velocities and Reynolds Stress 

* Growth Rates 

* Schlieren 
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- Goebel, S.G. and Dutton, J.C., “Experimental Study of Compressible Turbulent Mixing Layers,” AIAA Journal, vol. 29, 
no. 4, pp. 538-546, April, 1991. 

- Goebel, S.G. “An Experimental Investigation of Compressible Turbulent Mixing Layers,” Ph.D. Thesis, Dept. of Mech. 
and Ind. Eng., Univ. of Illinois., Urbana, Ill., 1990. 

-Gruber, M.R. and Dutton, J.C., “Three-Dimensional Velocity Measurements in a Turbulent Compressible Mixing Layer,” 
AIAA Paper 92-3544, July 1992 


Test Cases 
Compressible Mixing Layer 
Mean Velocity 
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Test Cases 
Compressible Mixing Layer 
Streamwise Turbulence Intensity 
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Test Cases 
Compressible Mixing Layer 
Transverse Turbulence Intensity 
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Compressible Mixing Layer 
Turbulent Shear Stress 


Test Cases 
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Test Cases 
Compressible Mixing Layer 
Shear Layer Thickness 


Shear layer thickness definition: 
The distance, b, between transverse locations where: 


ii = i, — 0.1A@ and ii = i, + 0.1A%. 
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Test Cases 
Compressible Mixing Layer 
Summary of Results 


¢ Results using FUNS3D with the SST-V model agree well with the Wind-US SST-V results. 
¢ All of the models compute the velocity profiles in the mixing layer well. 


¢ The Wilcox and SSG/LRR RSM and the EASM turbulence models are better than the SST-V 
model at predicting the turbulence quantities u’u’, v’v’ and u’v’. 


¢ The Wilcox and SSG/LRR RSM models give very similar results for v’v’ and u’v’. For u’u’, 
the Wilcox RSM model does slightly better. 


SIGNIFICANCE 
The Wilcox and the SSG/LRR full Reynolds stress turbulence models give improved turbulence 
predictions over the SST-V two equation turbulence model for this supersonic mixing layer case. 
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Test Cases 
Axisymmetric Subsonic Jet 
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- Bridges, J. and Wernet, M. P., "Establishing Consensus Turbulence Statistics for Hot Subsonic Jets," 
AIAA Paper 2010-3751, 16th AIAA/CEAS Aeroacoustics Conference, Stockholm, Sweden, June 2010. 
- Bridges, J. and Wernet, M. P., "The NASA Subsonic Jet Particle Image Velocimetry (PIV) Dataset," 


NASA/TM-2011-216807, November 2011. 
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Test Cases 
Axisymmetric Subsonic Jet 
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Test Cases 
Axisymmetric Subsonic Jet 
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Test Cases 
Axisymmetric Subsonic Jet 
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Test Cases 
Axisymmetric Subsonic Jet 
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Test Cases 


Axisymmetric Subsonic Jet 
Summary of Results 


¢ The SSG/LRR model shows some benefits over the SA and 
SST-V models at predicting the mixing. 
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Summary NASA 


Two RSMs available in the FUN3D code, the Wilcox and the SSG/LRR, were 
evaluated for four test cases: a transonic diffuser, a supersonic axisymmetric 
compression corner, a supersonic compressible planar mixing layer, and a subsonic 
axisymmetric jet. 


RSM results were compared with solutions computed using the SA and SST-V 
turbulence models, and an EASM (planar mixing layer). 


Transonic diffuser - results were somewhat inconclusive as to the benefits of the 
RSMs 


The supersonic axisymmetric compression corner — the SA model was best for 
computing the pressure rise and the separation location and length. The Wilcox RSM 
gave results similar to SST-V, and the SSG/LRR RSM severely over-predicted the 
onset of separation. All models had difficulty computing the boundary layer profiles 
and turbulence quantities in the separated region, and no additional benefit was 
gained by using RSMs. 


Supersonic planar mixing layer — the RSMs gave the best predictions of the 
turbulence intensity, turbulent shear stress and shear layer thickness 


Subsonic axisymmetric jet - SSG/LRR predicted the mixing of the core velocity the 
best 
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Conclusions 


The four cases examined are flows that are challenging for current 
turbulence models because they contain mixing, shock waves and/or 
separation. 


Overall, the RSMs showed benefit over the SA and SST-V models for the 
planar mixing layer and the axisymmetric jet flow, and may be useful for 
future nozzle calculations. 


While the cases examined are challenging flows, they are still relatively 
simple in geometry and flow features. 


More complex flow cases may reveal more benefits of the RSMs and are 
recommended for future study. 
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Questions? 


Extra Slides 
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Turbulence Models, cont’d 
e Seven-equation omega-based full Reynolds Stress models: 
¢ SSG/LRR RSM: 
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Turbulence Models, cont’d 
SSG/LRR RSM 
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Turbulence Models, cont’d 
SSG/LRR RSM 


Blending equation for @ = a, Bu, Ou, Oa - 


=F, o%+1-F)o© F, = tanh(¢*) 
c Vik 500u to ak 
= min|max | ——, - = 
C,wd pwd? }’ a) p 0k Ow 


ees ee, 2 
i) max (se ax, ,0) d 


Blending and closure coefficients for SSG/LRR RSM 
| te | ba — ee a a 


eee Te PPE eee 


al i a Kal la 
(€) 


40 


Turbulence Models, cont’d 
Wilcox RSM 
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Turbulence Models, cont’d 
Wilcox RSM, cont’d 
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Test Cases 
Sajben Diffuser — Strong Shock Case 


Wind-US and FUNS3D Results 
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Test Cases 
30° Axisymmetric Compression Corner 


Mach Contours 
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